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ABSTRACT 


Amagnetohydrodynamic flow of incompressible viscous fluid past a steadily moving infinite vertical porous plate is discussed under the action of transverse magnetic 
field. Considering a heat flux within the medium due to the constant temperature gradient at the plate, the effects of induced field and fluid concentration on flow 
and temperature are discussed for smaller values of Eckert number. The expressions and analytical solutions for fluid velocity, temperature, fluid concentration and 
induced field are obtained following Raptis [1]. The results are shown graphically followed by discussion and conclusion. 
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1. INTRODUCTION 

Problems of free and forced convective flow under the action of induced magnetic field are frequently encountered in various physical and engineering fields. The 
phenomenon of natural convective flow is not often caused entirely by the effect of temperature gradients but also due to the applied magnetic field which in turn induces 
another magnetic field that influence the flow largely. The study of magneto-hydrodynamic flow in the presence of porous walls has many applications in various fields 
of science and technology. Soundalgekar& Bhatt [1] have studied the oscillatory MHD channel flow and heat transfer. Soundalgekar & Bhatt[2] have studied the 
oscillatory MHD channel flow and heat transfer. Raptis & Soundalgekar [3] have studied the MHD flow past a steadily moving infinite vertical porous plate with mass 
transfer and constant heat flux. Drake [4], Soundelgekar [5], Pop [6], Dube [7], Singh [8], Das and Sanyal [9], Mana [10], Seth and Banerjee [11], Mittal et al.[12], and 
many other authors have studied various problems of MHD flow taking induced magnetic field into account. Singh and Singh [13] has studied on MHD effects on flow 
of viscous fluid with induced magnetic field. 


Motivated by the above referenced works and the numerous possible applications in the field Science & Technology specially in Geological Science, Astronomical 
Science and Biophysical Science, it is our interest to investigate a steady, fully developed MHD convective heat and mass transfer problem of an incompressible fluid 
flow where an incompressible electrically conducting viscous fluid moves past a steadily moving infinite vertical porous plate under the action of uniform magnetic 
field that induces another field transverse to it. The expressions and analytical solutions for fluid velocity, temperature, concentration and induced field are obtained 
following Raptis [1]. The effects of flow parameters like Schmidt Number, Hartmann Number, Prandtl Number, Magnetic Prandtl Number on fluid velocity, temperature, 
magnetic field, mass transfer, skin-friction and rate of heat transfer, are analyzed graphically and compared to those in absence of induced field. The results observed 
some significant differences in flow pattern in variation of flow parameters. 


2. FORMULATION OF THE PROBLEM 

The flow of an electrically conducting, incompressible viscous fluid is considered. The plate which is vertical and porous is assumed to be moving steadily in the vertical 
up- ward direction along which the x-axis is chosen. The y-axis is chosen normal to the plate. Under these conditions, the basic equations of combined free and forced 
convective flow,under usual Bossinesq’s approximation are derived by Raptis & Soundalgekar [3]. 


3 ASSUMPTIONS 


(i) The level of concentration of foreign masses are very low such that the Soret and Dufour are negligible. 

(ii) The flow is a free and forced convective flow under Bossinesq’s approximatnion. 

(iii) All the variables (i.e. velocity field, magnetic field, temperature field and mass concentration) are only dependent y-axis. 
(iv) Hall Effect, Polarization effect, Buoyancy effect and the effects of heat radiation are negligible. 

4. MATHEMATICAL FORMULATION 


We consider a fully developed laminar unidirectional flow of an incompressible, electrically conducting incompressible fluid past a vertical porous plate. The plate is 
assumed to be moving steadily in the infinite vertical upward direction along which x-axis is chosen, while y-axis is normal to the plate. Let a uniform magnetic field 


Hy) is applied along the motion of the plate that induced magnetic field |B =. H(y) transverse to the motion of the plate. 


5. GOVERNING EQUATIONS 
Under above assumptions, the velocity and magnetic field magnetic field components are 


V=[u(y), v0.0]; B=[H(y), wH(y),0 ] () 
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H &u, are being magnetic flied applied and magnetic permeability of the medium. 


Under above assumptions and the geometry of the problem, the governing equations are as follows. 
Equation of conservation of momentum 


Voge = v(@u/ dy’) + uH(y Sr + EB (T-T.)+¢ B(C-C.)} Q) 


Equation of thermal diffusivity 


pCp (Vor DEKT/ dy) +u CO) 2 Ey 3) 
Equation of Mass diffusivity 

q@ £)=De&> (4) 
Equation of Magnetic diffusivity 
@=VX(vXB)+ *VH (6) 


Introducing non-dimensional quantities as given below 


w=> ‘ y' = Vo- C’=(C-C.) /(Cw-C.); T’ = (CT - T. )k Vol / (qv) ; 


Vo 
H’= (H/v) V(pe / p ), qis a quantity with dimension M T + (6) 


Substituting above non-dimensional quantities and removing primes, equations (2) — (5) are as follows. 





2 

- ~ 6,0+6,7+24mM4 (7) 

ay dy? ay 
1, 07H OH au 
mal Oy? }+ oY +M,=0 (8) 
ac. 1 7d7C 
ar ke. lage (9) 
eT yar du \7, 1 (H\> 
ateg = PE aay) (10) 
where, 


v= * Kinematic viscosity ; = Coefficient of viscosity of the medium ; T = Fluid temperature ; 8 = Volumetric Coefficient thermal expansion ; B = Volumetric 


Coefficient of mass transfer ; « = thermal diffusivity ; Um = Heo (magnetic dif fusitivity) ; o=electrical conductivity, |. magnetic permeability of the 
medium ; C = fluid Concentration, B = Volumetric Coefficient thermal expansion ; f$ = Volumetric Coefficient of mass transfer ; D = Coefficient of mass transfer ; S. 
= Schmidt Number; M= (H/vo) V(ue/p), Magnetic Hartmann Number, 


Pi, = Heov Magnetic Prandtl Number; G, = (g B q 0?) / (k vo*), Grashoff Number ; 
Gm= [g Bv (Cw — Cx)/ (vo3), Modified Grashoff Number; 
P, = (uC, / k), Prandtl Number; E = (k vo*)/(quC,), Eckert Number 


The corresponding non-dimensional boundary conditions are 





aT 
y= 0, u=U, Gath C=1, H=K 
& yoo, ur0, T>0, C >0, H >0 (11) 


Where U and K are two assumed values for u and H respectively suitable to the problem 

6. SOLUTIONS OF GOVERNING EQUATIONS 

Solution for Mass Transfer equation under above boundary conditions is 

c(y) =e" Se» (12) 


In order to get solution of the governing equations (7 — 10) we consider following series expansion in terms of E, the Eckert Number. 


u=up + Eu, +0(’) 
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H= Hy) +E H, + 0(£2) 
T=T) +ET, + 0(E’) (13) 





Using above expansions the governing equations can be divided in increasing order of E (i.e. 1°, 2"! ----) as follows. 


1“ order of equations 

















07 ud uo OHo = 

Se = So MS + GTO HG C= 0 (14) 
0?T» OT _ 

(1/P) “P+ 2B =0 (15) 
0*Ho OHo iS 

(1 /Pm) apt 5 (1+ M)=0 (16) 


Corresponding non-dimensional boundary conditions are 

















y= 0, uo =U, @ )=-1, C=, Ho =K 

& you, u-—d0, T>0, C 750, H 50 (17) 
2" order of equations 
Oru, Ou, On, = 
mr “a tM a +G,T;+GnzC=0 (18) 
er; OT, duo \ 2 1, 0Ho \> 
pe ag gg a ay (19) 

07H, OH, du, 
(1 /Pm) ar ee =0 (20) 
7. VISCOUS-DRAG AND HEAT-TRANSFER 


The physicalquantities of ourinterestat the plate are Skin friction coefficients [t]y-» and rate of heattransfer [Q]y-9 areas follows. 


1, aT 


t=-u(),<0 and Nu =a ei ly= 0 


Using (6), the non-dimensionalviscous drag per unit length i.e. Skin friction and rate of heattransferat the plate y = 0 are given as 


Tp =(F)s-0 (21) 
Ni =15 Iso (22) 
Where, 


d aT ‘ : : F rey ’ 
(=) and (=),are the non-dimensional gradients of velocity and temperature with in the medium. 


8. SOLUTION OF GOVERNING EQUATIONS. 
The solutions of governing equations (14)to (20)under their boundary conditions are as follows. 

me atta Pp 
To= me y (21) 
Hy =Cye™+ (ks_c,)e™” —kze?m¥+ kyeSc¥ (22) 
Up=aye™Y + a e™Y — azeS¥+ aye? mY (23) 


ky = Pat, ka = Pm(M? — 1), 

















‘ MP nc, 

Pm (Pm? sa kiPm a k2Pm) 

ee PmMGm 
ss —S5.3 — ky, + kaS¢ 
_ m2 = mM, = = 
ke = (we 1) (me 4: 1), kg =K—(ky ks) 
m Ss P 
key = kes (7 +1) +k, (145% )+i (14+ ) 
m m m 
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ky + fia? + 4k, —ket ka? + 4k, 





Ms,6 = 
























































. 2 : 2 
% = 2 
Kg=p,+1 > Ko= P,(1—M ) 
ae —kg + Vkg? — 4ko 
3,4 2 
Utky — 7m Utky 
a= ke C2=ks — Cy ay = Gas 1)¢ ke ); 
U+k7 M2 Sc Pr 
dy = (ks — 9) GE 41), ag = ka (1-) . ae = ea + PE 
6 m Pm MPm 
ey Eye Sa iD Da Go ae ai yaar Dey = — _ = 
Ag = AM{; Ag = AZMZ 3 A7 = 2A,AZM,Mz 5 Ag = A3S¢Ag = AZP{, Ayo = 2AgQ4S-F., Ay, = 24,a3M,S_,ay2 = 2aya,M,P,,a,3 = 12aza3mz, 
_ = —P,as _ —Prdg _ —Praz _ —Prag _ —Prdg _— _ —Praio 
ayg = —2a2a4MzPy, ars am2—P,2m,"916 — am2—P,2m, 917 ~ Gm2—P,2m, 918 ~ Gm2—P,2m,? 219 ~ am2—P,2m,? 220 ~ Gm2—P,2m,’ 
si Tr a 1 Tr 1 1 Tr 1 1 © =f 1 Tr 3 x i a 4 
a = —Pray1 a a, —Pray2 a au —Prai3 a = —PrAy4 a pia? MPmG, a a, MPG, 
21 ~~ 4m2—P,2m,’ “22 ~ 4m?=P,2m,’ “23 ~ 4m?2—P,.2m,’ “24 ~ 4m?=P,2m,’ ©25 ~ (2m,)3+ke(2m,)2+k7(2m,)’ ©25 ~ (2mz)3 +k g(2mMz)2+kz(2Mz)’ 
MP mGr MP Gr MP m Gr MP m Gr 
O97 a et a DS a Poe a2: DAS DAZ 80 eT DO 
(™14m2) +ke(m4+m2)*+k7(m, +m) (-2S¢)? +k6(-2S¢)?+k7(-25¢) (2P,)° +k (2Pr)* +k7(2Pr) (Pr-Sc)? +6 (Pr-Sc)* +7 (Pr-Sc) 
MP m Gr MP mGr —MPmGy 
O34 ena ia OD Bk es LAR ED? SS ee EONS Le Den 
(m1 -Sc)? +6 (M4-Sc)* +k7(M4-Se) (m,+P,)? +ke(m4+P-)*? +k7 (m1 +P,) (mM2-Sc)? +K6(M2—-Sc)* +k7(M2-Se) 
—MP Gy —P,.b? 
az4 = by = CyM3, bz = (Ks_c,)m,, bz = —k3P., by = —k,S,, b¢ = ——S"—~ 
34 = Gay eP Ptkglng+P,)#R CmgaP,» Pt = CAMs, By = (Ks—c)imy, Bg cae aie 40? 86 Pm (43 —P,2M3)’ 
pe —P,b2 _ —P;-2b,b2 
ig Pm(4m3—P;-2m4) bs Pm{(m3+m,)?—P;-(m3+m4)} 
ss —P,b2 _ =Prbi 
bo = Pi@P,)2—P.2P)) ho =O 
m{(2Pr) —P,(2P,)} Pm[(—2S¢) —P,(—25S¢)] 
ae —P,2b3b4 a —P,2b,b3 — —P,.2b,b4 eo —P,2Db3 b= —P,2D,D4 
11 pL (Py—-Sc)2—(Pp—-Se) Pr]? 12 Pr [mg +P,)2—P,- (3 +P,)] 23 Px [Cm3—Se)2—P-(m3—-Se)] 84 Pl (my +P-)2—P-(m,tPy)]) 22 Pn [(m4—-S_)2—P, (M4 Se) 
m T c T c. Tr. m 3 Tr. Tr 3 Fr m 3 C. r 3 C. m 4 Tr. Tr 4 F m 4 (32 r 4 S 
Pie MP Gy i MP Gy ice MP Gy Peek MP Gy ies 
16 ~~ (2m4)2+K¢(2m4)2+Ky (2 ue eLh 24K k > 818  (2P,)34+k6(2P,)+k7(2P,) > 19 (—25,)3 +h (—25,)+kz(-2S,) > 29 
(2m,)?+k6(2m,)* +Kk7(2mM,4) (m3t+m,4)? +k_(mz+m4)+k7(m3+m,) (2P,)3 +k 6 (2P,)+k7(2P,) (-25¢)? +k¢(—25S-)+Kk7(-25¢) 

MP Gr = MP Gr = MP Gr a MP Gr - 
oe bn = oo ba = oo ba DoD by = 
(Pr-Sc)? +6 (Pr-Se)+k7(Pr-Se) (m3+Pr)? +ke (m3 +P,)+k7 (m3 +P) (m3-Sc)*? +K6(mM3-Sc)+k7(M3-Se) (m3 +P)? +Kkg(m3+Py)+k7(m3 +P) 

MP Gr 
(m3-S¢)3 +K6(m3-Sc)+k7(m3-Se) 

T, = cgePr? 4 ayge™ 4 ayge™” 4 ayze(™*™DY 4 ayge 25" 4 aygePr? + a,, grsoy + az1e(™ SY 4 agreO™*PY _— qyge(M2-Sc)Y _ 


Ape 24 Pr)Y he e2M2Y be e2M4Y he (MstMa)Y , be e2Pr , bye ~25cY yyy ge Pr-Sc)Y , by e M3 t PHY by eM SV ,h,ge(Mat PY ,h, ,e(Ms-SO)Y 


Hy=cge™®” +c7e™ scgePr” sap 5e7™Y 40260729 4Ap7¢ mr+ma)y +028Azge °° + Az9e2P"Y ,agge Pr-Sc)Y 434 COM “SCY g,e MIMAY Gz,e(M2-Sc)Y — 


340024 PLY by e2™3Y by OY spe M3tMYY by,  e2PrY sig 4e 25) gDo ge PrS)Y by e M2 tPr)Y bo e (M3-S)Y jy9393e(M3tPr)Y -ho,e(M3s-Sc)Y_ 


1 es 

uy = — hipy (Comse™™” + c7mee™e + CgprePrY + Arg2mye?™Y 4 azg2my2e7™ 4 az7 (My +Mz eOMA*MY 4 gg (-28c)@ 78 4 Agg2e7PrY 4.39 (Pr - 
m 

Seri) e(Pr-Sc)¥ + Qa, ( Mm - Se ) e(mi-Se)y + A327 ( Mm + Pr ) e(™itpr)y - 33 ( My, - Se) elm Sc) 


G34 (Mz +P, ) e™2tPY . bye2mze2™3Y 4 bygmye2™4 4 by7 (Mz +My) COM3tMY , by g2p,e2PrY 4 Dyg (25, C075) 4 Dog (Py - Se) CPPSAIY 4 
bay (Mg + Py) eM*PHY 5 bop (Mg - S_ EO S)Y 5 bog (Mg + Py COM *PHY 4 bog (Mg - Sc CM -SO)Y -= | Cge™5Y + c7e™6Y 4 CgePY 4 a25e™ 
+ Agge?™Y 4 AyzeMtM)Y , arge~*c¥ 4 Agge2PrY 4 AggePr-SI)Y¥ 4. agyeO™~S)Y , agye(™+PrY - Qgze(MSd)Y - qz,e(M2tPr)Y 
+b 5e2™2Y ,b, .e2™4Y by 7e Mt MY bh, ge2PrY by ge 25) sb ge Pr-S0)¥ .hy,e(M3t Pr)Y ho e(M3-SC)Y , by ge M3tPr)Y bye (Ms—Se)Y) 


((Cg + Qg5 + Az6 + Az7 + Agg + Azgg + Azo + Ag, + Azz - 433 - Aga + dys + Dye + yz + big + dio + bao + bg, + 


C741 
baz + bo3 + b24)) Ms 

C72 = CePy + Ag5 2 My + Azq 2 Mz + Ag7 (My+ Mz)+ Agg (-2 Se): 

C73 = 2 PrA29 + A390 (Pr - Sc) + Az1 (My, - Sc) + G32 (My, + Pr) - A33° (M2 - Sc)- Azq (Mz +Pr)+ dys 2 mM3t dye 2 Mg + dy7 


(m3 + M4) + big 2 Dr: 
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C74 = Dbyg (-2 Sc) + bag (Pr - Sc) + bar (M3 + Pr) + boa (M3 - Sc) + b23 (M3 + Pr) +b24 boq (M3 - S¢): 
C7 =(1/ (ms- mMe)) (C72 + C73 + C74 - C71): 
C6 = -(C7 +c8 + A25 + Azg + Az7 + Azg + Azgg + Azq + Azz + Az2 - AZ3 - Aza + bys + bi6 + by7 + big + by9 + boot bo4 + b22 + bo3 + 


b24)) 


9. RESULT AND DISCUSSION 
The purpose of study is to discuss the effects of induced field (p,.) and mass transfer (S,) on velocity field(u), temperature field (T), magnetic field (H), skin friction at the 
plate [t ]y-o, and the rate of heattransferat the plate [Q],-0 . In calculating numerical values of the physical quantities u, H, T, a 7 oe tand Q, we have considered 


Gm = 4 (arbitarily); E=0.001 so to satisfy the series expansion [13]; Gr = 10 because it relates to the various problems in engineering and technologycal activities; 
Mass transfer parameter Sc varied as 0.22 (for Hydrogen), 0.6 (for Water vapour ) & 0.75 (for Oxygen). The physical parameters whose effects on flow are the 
objectives of this study i.e the induced field parameter Pm (Magnetic Prandtl Number )is varied as 0.5 to 4.5 while applied field parameter M (Hartmann Number) is 
varied from 5.0 to 20.0; We consider that Pm < 0.5 means slow affect of induced field while Pm >= 2.5 moderate affect of the induced field. Further the parameter 
relates viscous with thermal diffussion i.e. Prandtl! Number (Pr) is varied as 0.015 ( for mercury) , 0.16 ( for noble gase ), 0.7 ( for water at 20° c ).Following are the 
results from the figures (1-5) drawn using numerical values of the physical parameters given above. 





Fig. 1(i— ix): Velocity Profile, variation of u with Pm for Sc = 0.22; Gr=10; Gm =4; U=1; K=.5; E=0.001; 
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Fig 1. Shows variation of fluid velocity with Pm, M & Pr. It is clear from the figures that when M remains unchanged, u increases with Pm; the rate of increase is higher 
for higher values of Pm. keeping Pm constant, and M rises up, u increases with largely; e.g. at y = 1.0, & Pm = 1.5, u = 815.7 for M=5.0 while u = 22581.4 for M = 10.0; 
similarly for Pm = 0.5 and Pm = 2.5. For higher values of M (> 10) and lower value of Pr, velocityu reverses its direction. The variation of u with Pm almost same in 
nature when Pr increases, keeping M unchanged. 
























































Fig 2 (i— ix): Temperature Profile, variation of T with Pm for Sc = 0.22; Gr =10; Gm =4; U=1; K=.5; E=0.001; 
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Fig 2. Show variation of fluid Temperature (T) with Pm, M & Pr. The figures show that within the medium, T decreases with the rise of Pm for all values of M (= 5.0, 
10.0, 15.0) & Pr (= 0.015, 0.16, 0.7). Keeping Pm constant, T decreases with increase of M largely(e.g. at y=0.8& Pm=1.5, T=-100.26 for M=5.0 while T = - 


3968069.0 for M = 10.0; similarly for Pm = 0.5 and Pm = 2.5 ); this implies higher rate of heat transfer el at higher values of M. When M & Pm are at constant, T 


decreases with increase of Pr ; this implies lower the thermal diffusion, lower fluid temperature. The nature of variation of T with Pm atconstantM is almost same 
for different higher values of Pr (= 0.015, 0.16, 0.7). 











Fig 3 (i — iii): Magnetic field Profile variation of H with Pm, for Sc = 0.22; Gr=10; Gm =4; U=1; K=.5; E=0.001; 
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Fig 3. Show variation of magnetic field with Pm, M & Pr. Away from the plate, H decreases for all values of M,Pm & Pr ; the rate of decrease is more when Pm is 
higher; similarly for M when Pm is constnt. Field a very near to the plate is less variant of Pm depending upon applied field (e.g. field is less variant within y<2.0 
when M = 5.0 while itis y <5.0 when M= 10). When M increases at constant Pm, H decreases largely (e.g. at y=0.4,& Pm=1.5; H = -21.022 for M=5.0 while H 
= -162.344 for M = 10.0 and H = --1437.68 for M = 15.0); similarly for Pm = 0.5 and Pm = 2.5. 
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Fig 4 (i — ix): Shear Stress (t ), variation with Pm & M_ for Sc =0.22; Gr =10; Gm =4; U=1; K=.5; E=0.001; 
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Fig 4. Show variation of Shear stress (t) experienced by the fluid motion at the plate ( y=0 ) with Pm, M & Pr. t at the plate y = 0, varies inversely with Sc when M, 
Pm &Pr are fixed. For all values of M, t decreases rapidly with the increase of Pm within lower values of Pm <2.0 (approx.), while increases slowly beyond it. When 
Pm, M &Pr are constant, t increases with the rise of M ; decreases with rise of Sc at constant Pm, Pr& M. When field is high M =20, the variation of t with Sc 
becomes less significant but becomes prominent when both M &Pr are high. This implies the strong effect of Pron t. 
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Fig 5: Rate of heat transfer, proportional to ( 5 ) Variation with Pm & M for Gr=10; Gm =4; U=1; K=.5; E=0.001; 
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Fig 5. Show variation of Q away from the plate. A little away from the plate y >O (depending upon M), Q is less variant for the variation of Pm and Pr (e.g. for M= 5.0, 
y < 4.0 aprox. while M= 10.0, y < 2.0 aprox.); beyond this, Q increases sharply with the increase of Pm. The limiting value of y beyond which Q varies linearly with Pm 
depends upon M. The rate of increase of Q within the medium is more when Pm is higher. The variation of Q within the medium is higher when either of M &Pr is 
higher and when both M &Pr are higher Q is much higher within the medium. 


10. 


11. 


CONCLUSION 


Applied field(H), Induced field (Pm) and thermal diffusivity factor ( Pr) have influence over fluid velocity. Fluid velocity increases with the increase of applied 
field and induced field. The direction of fluid velocity reverses when M & Pm is higher but Pr_ is lower. 

The effect of H, Pm ,Pr on Fluid temperature (T) is opposite in nature to the variation of velocity as above. Temperature is much low when all the three 
parameters Pm, M &Pr are higher. This implies strong thermal diffusivity effect with magnetic interaction and also higher rate of heat transfer within the 
medium. 

Away from the plate, Magnetic field within the medium decreases with the rise of induced field; this implies strong magnetic diffusivity within the medium 
away from the plate. At very close to the plate, magneric field is almost uniform (i.e. less variant with induced field) depending upon applied field. 

The effect of applied field on shear stress significantly depends upon induced field. The nature of variation of shear stress is almost opposite within lower 
level of induced field and beyond it; the level depends upon applied field. Bothinduced field and mass transfer parameter have significant role over the shear 
stress. 

Applied field, Induced field and thermal diffusivity factor have distinctive influence over the heat transfer. Heat transfer is more at higher induced field. It is 
further higherwhen thermal diffusivity factor (Pr) is higher along with higher applied field (H)but less induced field(Pm); this implies higher rate thermal 
and magnetic interaction. 
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